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Abstract 

The perturbation approach is used to derive the exact correlation length 
£ of the dilute Ax lattice models in regimes 1 and 2 for L odd. In regime 2 
the A3 model is the Eg lattice realisation of the two-dimensional Ising model 
in a magnetic field h at T = T c . When combined with the singular part f s 
of the free energy the result for the A3 model gives the universal amplitude 
f a £ 2 = 0.061 728 . . . as h -> in precise agreement with the result obtained by 
Delfino and Mussardo via the form-factor bootstrap approach. 



The integrable E 8 quantum field theory of Zamolodchikov ffl, H is known to be in 
the same universality class as the two-dimensional Ising model in a magnetic field at 
T = T c . Moreover, an integrable lattice realisation of the Eg Ising model is provided 
by the dilute A3 model || |4[], upon which explicit exact and numerical calculations 
pertaining to the Ising model in a magnetic field can be performed [3-13]. 

In this letter we present the correlation length of the dilute Ax lattice models in 
regimes 1 and 2 for L odd, for which the off-critical perturbation is magnetic-like. 
This includes the magnetic correlation length for L = 3, of relevance to the magnetic 
Ising model at T = T c . 

The dilute Ax model is an exactly solvable, restricted solid-on-solid model defined 
on the square lattice. Each site of the lattice can take one of L possible (height) values, 
subject to the restriction that neighbouring sites of the lattice either have the same 
height, or differ by ±1. The Boltzmann weights of the allowed height configurations 
of an elementary face of the lattice are |3|, 

#i(6A-u)tfi(3A + u) 




tfi(6A)tfi(3A) 
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The crossing factors 5(a) are defined by 
and ^i('u), i?4(u) are standard elliptic theta functions of nome p 

oo 

0i (u) = 7?i(w,p) = 2p 1/4 sin« [] (l -2p 2n cos2 M ip 4n ) (l - p 2n ) (3) 

n=l 

4 (u) = ^4 («, p) = ft ( X - V n_1 cos 2« i p 4 "- 2 ) (l - p 2n ) . (4) 

n=l 

In the above weights the variable A and the range of the spectral parameter u are 
given by < u < 3 A with 

A = - 7T (5) 

r 



where r = 4(L + 1) and s = L in regime 1 and s = L + 2 in regime 2.[] The magnetic 
Ising point occurs in regime 2 with A = 

The row transfer matrix of the dilute A models is defined on a periodic strip of 
width N as 



T {b} 

1 {a} 



N 



bj b j+1 
a,j dj+i 



where {a} is an admissible path of heights and cln+i = cti, bu+i 
we take N even. 

The eigenvalues of the transfer matrix are || [14], |15| 
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where the N roots Uj are given by the Bethe equations 
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and uj = exp( i 7r£/ (L + 1)) for £ = 1, ... , L. 

There are several methods at hand to calculate the correlation length. Here we 
apply the perturbative approach initiated by Baxter [IB, [T7[|. For L odd this involves 
perturbing away from the strong magnetic field limit at p = 1. We thus introduce 
the variables 

w = e - 27ru/e and x = e~ 7r2/re (9) 
conjugate to the nome p = e~ e . The relevant conjugate modulus transformations are 
'7r\ 1 / 2 „ (1Q) 

(11) 
(12) 
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where q = e _7r /e and 
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E(z,p) = \[(l-p n - l z)(l-p n z- 1 )(l-p n ). 
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In the ordered limit (p — ► 1 with w/e fixed) the Boltzmann weights for L odd 
reduce to 



d c 
a b 



w H( d ,a, b) , 



(13) 



1 The model has other regimes, but they are not of interest here. 



The function H(d,a,b) is given explicity in ||, being required for the calculation 
of the local height probabilities. In this limit the row transfer matrix eigenspectra 
breaks up into a number of bands labelled by integer powers of w. In regime 1 there 
are |(L + 1) ground states and in regime 2 there are \{L — 1) ground states, each 
with eigenvalue Ao = 1. The bands of excitations are relevant to the calculation of 
the correlation length. 

The number of states in the w band is \{L — 1)N in regime 1 and \{L — 3)N in 
regime 2. These correspond to introducing in all but one of the ground state paths 
{a} a single non-ground state height, in any position. In particular, note that there 
are no excitations in the w band for L = 3 in regime 2. Thus for the magnetic 
Ising model we must consider excitations in the w 2 band. These are harder to count, 
arising from a variety of both single and multiple deviations from ground state paths. 
However, we observe numerically that (apart from when N = 2) there are 4N states 
in the w 2 band. 

We associate a given value of £ with each eigenvalue by numerically comparing 
the eigenspectrum at criticality (p = 0) with the eigenspectrum of the corresponding 
O(n) loop model |1| for finite N.f\ Each eigenvalue can then be tracked to the ordered 
limit. The band of largest eigenvalues is seen to have the values £ — 1, . . . , |(L + 1) 
in regime 1 and £ = 1, . . . , \{L — 1) in regime 2. 

Setting Wj = e~ 27TUj ^ e , the eigenvalues (|7|) can be written 
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The Bethe equations 

E(x 2s /wj,x 2r ) 
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The calculation of the largest eigenvalue proceeds from the x — > limit with w 
fixed in a similar manner to that for the eight-vertex [|16| and CSOS (T7| models. 
Each of the degenerate ground states has a different root distribution {wj} on the 
unit circle, depending on £. Defining the free energy per site as / = iV _1 log A our 
final result is 



cosh[(5A — 7c)7ck/e] cosh(irXk/e) sinh(7ru/c / e) sinh[(3A — u)%k/e] 



(16) 



. I k sinh(7r 2 /c/e) cosh(37rAfc/e) 

in agreement with the previous calculations via the inversion relation method [§, d, f| • 



2 Strictly speaking we compare with the eigenspectrum of the corresponding vertex model with 



In regime 1, the leading eigenvalue in the w band has £=~(L + 1) + 1. The 
root distribution has N — 1 roots on the unit circle and a 1-string excitation located 
exactly at wn = —x r . Applying perturbative arguments along the lines of [T7|] yields 
the leading excitation in the w band to be 
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At the isotropic point w = x 3s this reduces to 
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For L = 3 in regime 2 extensive numerical investigations of the Bethe equations 
have lead to a convincing conjecture for the thermodynamically significant strings 
@, 0. We find that the leading excitation in the w 2 band is a 2-string with 1 = 2. 
However, the state is originally a 1-string for small p. Such behaviour has been 
discussed in ||. Tracking this state with increasing p reveals that the 2-string is 
exactly located at — x 1 in the limit p = 1. There are finite-size deviations away 
from this position for small N and < p < 1. The location we find for this string is in 
accord with the previous numerical work |6|, . Applying the perturbation arguments 
in this case yields the leading excitation in the w 2 band for L = 3 to be 



A 2 _ 2 E(-x/w, x 60 ) E(-x n /w } x 60 ) E(-x 31 w, x m ) E(~x 41 w, x 60 ) 
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At the isotropic point w = x 15 this reduces to 
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The correlation length £ can be obtained either by integrating over the relevant 
band of eigenvalues or via the leading eigenvalue in the band at the isotropic point 
(see, e.g., P7| ). Doing the latter we have 



-i 



, A 
log-r- 



where A is the relevant leading eigenvalue. Our final results are thus 
for L odd in regime 1, with 
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for L = 3 in regime 2. 

The derivation of the correlation length for L ^ 3 in regime 2 is complicated. 
In this regime the leading excitation in the w band has i = |(L — 1) + 1 and, like 



the leading 2-string in the w 2 band for L = 3, it begins life for small N and p ~ 
as a 1-string. We have not pursued this further. Nevertheless we have numerically 
observed that the final result ( |TTD also applies to the leading w band excitation in 
regime 2. We thus believe that the correlation length fl22|) and the corresponding 
exponents below also hold in regime 2 for L ^ 3. 

It follows from (22) that the correlation length diverges at criticality as 

Z~JLp-"» as p^O (24) 

where the correlation length exponent Uh is given by 

!2(L+1) ■ i 

3L re S lme 1 

(25) 
2(L+1) • o 

3(1+2) re S lme 2 • 

The correlation length exponents are seen to satisfy the general scaling relation 
2vh = 1 + 1/5, which follows from the relation 

f s £ 2 ~ constant (26) 

where f s ~ p 1 + 1 /< 5 is the singular part of the bulk free energy and the exponents S are 
those following from the singular behaviour of ((T6|) || f|, |]| 

57T 

1(5 



The magnetic Ising case at A = ^ is of particular interest. From (P]) we find 



sin 

/ s ~4a/3 5 -p 16/15 as p^O. (27) 

cos 30 

On the other hand, from (|23| ) we have 

t - a J ■ , P' 8/15 ^ p^O. (28) 
8v3 sin | 

Combining these results gives the universal magnetic Ising amplitude 

/ s £ 2 = = 0.061 728 589... as p^O. (29) 

16 v 3 sin | cos ^ 

This is in precise agreement with the field-theoretic result obtained recently by Delfino 
and Mussardo, starting from Zamolodchikov's S'-matrix and using the form-factor 
bootstrap approach |19|, p0|. Full details of our calculations will be given elsewhere. 
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3 The same correlation length exponents should hold for L even, for which the integrable pertur- 
bation is thermal-like. The scaling relation is now 2v t = 2 — a, where v t and a are as given in ( |25| ) 
and H [|, ||, respectively. In particular, ( |25| ) gives the Ising value v t — 1 for L = 2 in regime 1, as 
expected. 
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